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PREFACE. 



My object in publishing these Notes on Finite DiflFerences is 
to give, in a convenient fonn, a collection of those elementary 
propositions a knowledge of which is required of the students 
who present themselves for the Institute of Actuaries* First 
Examination. I have not attempted a complete account of even 
that small portion of the subject utilized in actuarial science, 
thinking it desirable to confine myself almost entirely to the 
methods of elementary algebra. The only instances in which a 
knowledge of mathematics beyond the Binomial Theorem is 
required of the reader are § 2, Chap. IV ; Example 28 ; and the 
Note to Example 27. These are marked with asterisks, thus *. 

The difference symbol has been denoted by S when the incre- 
ment of the independent variable is taken to be unity, and by A 
when this increment is not so restricted. 

I am indebted to Mr. T. B. Sfbaoue, President of the 
Institute of Actuaries, for some valuable suggestions, and regret 
that I have been able only in part to avail myself of them. 



A. W. S. 

2, King William Street, 

London, M,C, 

20 Feb, 1885. 




§ 1. Aa the following short account of some of the more 
elementary theorems of Finite Difierencea ia intended foi" the 
use of studeiita who have no acquaintance with the methods of 
analytical geometiy, it aeema deairahle to preface it hy a short 
explanation of the method of representing variable quantities by 
curved lines or diagrams. 

§ 2. Let us eonsider the quantity a:'. Its vahie, of course, 
depends on that of ar, and if any particular ^-alue be given to a.; 
the coiTesponding value of j^ may be determined by multiplica- 
tion; e.g., a a:=3, 0^=9. 

Any quantity such as a;" which depends on x in such a 
manner that for each value assigned to x it takes a determinate 
value is called a function of x. 

As another illustration of a function, consider, out of 100,000 
persona born alive, the number, ti, who would be living at age x 
on the assumption that these persons were subject to some 
definite law of mortality. Here for each value of x, the age 
attained, wc have a definite number alive at that age. The 
quantity U is therefore, according to our definition, a function 



Any function of x may be conveniently denoted by the 
symbol Wjr. With reference to «,, the quantity x, (to which we 
may assign any value we pleaae, the corresponding value of «a 
being then determined), is called the independent variable. With 
reference to x, the quantity Ux, (whose value we eonsider deter- 
mined when that of .r is given), is called the dependent variable. 

§ 3. Suppose we wish to examine the series of values which 
a function takes for a aeries of different values of the variable .v 
which it involves. A very valuable and powciful method of_ 



conducting nach an examination, and of investigating the pro- 
perties of any function, is that of representing it by a diagram 
as follows : — 

Let a straight line, which we will call OX, be drawn on a 
sheet of paper from a definite point 0, and let lines measured 
along OX represent, according to some given scale, values of x. 
For example, if wc take one inch to represent the unit of x, 
x=5 will be represented by measuring along OX a line, ON, 
5 inches long. 

In order to represent the value of the function for any value 
of w, draw through the end of the line which represents the value 
of a: a line perpendicular to X, and of a length correMpondiug 
on the same, or, it may be, a dif- ,. 

ferent scale to the magnitude of 
the function for this value of x. j 

For example, if the function is ar^, I 

aud the scale in each case one to \ 

the inch, for x=2, and therefore " -"' ■*' 

«'=4, a line ON 3 inches long must he measured along OX 
and through the eud of it, N, a line A' P drawn perpendicular to 
X and 4 inches long. 

In this way, when the algebraic expression for the function is 
given, any number of lines such us N P may be drawn, each 
representing the value of the function for that value of x which 
is represented by the line between and the foot of the per- 
pendicular. The summits of these perpendiculars will form a 
aeries of points lying on a curved line, and the process above 
described is that of mapping out the diagram or curved line 
correapondiug to the function. 

§ 4. To complete the representation of any function of x it 
is necessary to provide for negative values. This is done by pro- 
ducing X through and maiing the convention that negative 
values of x shall be represented by lines meaam-ed along the 
produced line, which we will call X', and negative values of the 
function by drawing the perpendiculars down instead of up. 

It will he a useful exercise for the student to rasp out the 
diagrams for two or three different functions. 
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Take, for cxamplcj the function i—x'. The diagram will be 

found to take 

the form alio wn 

by the dotted 

curved line 

given in Pig. 

{2). The curve 

is that known 

as a parabola. 

For the point 

P, X and the 

function are 

both positive ; 

for P' J! is positive and the function negative, for P" x is 

negative and the function positive, for P"' both are negative. 

Lines such as O iV measured along X or OX' are usually 
called abscissie, and the perpendiculars, such as P N, ordinates. 

§ 5. Su])po8e we are given, not the algebraic expression for. a 
function, but the curved line or diagram which foi-ma its complete 
representation. Then, in order to find its value for any given 
value of a', we must measure along OX & line, say ON, repre- 
senting X in magnitude, and draw through N a line perpendicular 
to N, so as to cut the curve. If P be the point of intersection 
of the pci-pendicular and the curve, then N P will represent the 
value of the function corresponding to the given value of sc. If 
this perpendicular, which may be drawn both up and down, meets 
the curve in more than one point, there will be more than one 
value of the function for the given value of x. If, on the other 
hand, it does not meet the curve at all, there will be no value of 
the function for the given value of x. 

As an illustration : The curve corresponding to the function 
which is the square root of x may be drawn on a sheet of paper 
by u very simple and mechanical contrivance. Supposing it so 
drawn, then we can ascertain the square roots of numbers by 
merely measuring lines on a diagram. 

The form of the curve in question is shown in Pig, {3), For 
each positive value of x there aie two perpendiculars, P N P'N, 
equal in length but on opposite sides of OX, corresponding to 




the fact that every positive number has twu squarG roots equi 

maguitude but opposite in sign. 

§ 6. It will be useful to regard the lines iV and P N, 
which are called 
respectively the 
ahscisaa and or- 
dinate of the 
point P, from a 
different point of 
view. Starting 
with a given 
line of reference 
^OZ',the posi- 
tion of any point on the diagram is determined when its abscissa 
and ordinate are given. The absciasu we have denoted by the 
letter *, the ordinate is usually denoted by the letter y. Thus the 
point whose abscissa is 3 and ordinate 10 — in other words, fot 
which x=3, i/=10 — is the point P, found by measuring along 
OX a line N=S and drawing through N a line JVP per- 
pendicular to JV and equal to 10. This point might, however, 
equally well have been found as follows : — Take another line 
of reference J" 1" drawn through O at right angles to X X', 
measure along Y a line OSf=lO and through M draw M P 
perpendicular to F 
and equal to 3. We 
see, in fact, that the 
abscissa (denoted by iv) 
is the distance of the 
point P from the line 
FOr, called the axis 
of y ; and the ordinate 
(denoted by y) is the 
distance of the point P 
from the line X X'l 
called the axis ot x. 

Looked at in this way there is no distinction in theory between 
the abscissa and ordbiate of a point. The two are called the 
ca-ordinates of the point. 
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Afl an example, consider the function which i 
Foi- this function we have the equation 



the square of 



which may be called the equation to tlie curve which represents 
the function under consideration, since it is the equation which 
connects the co-ordinates of every point on this curve. Giving to 
the independent variable x any values we please, and determining 
from the equation the corresponding values of y, we may obtain 
any number of points on the carve. We may, however, proceed 
differently. The equation may be WTitten 

and we may now take y for the independent variable, assign to it 
any values we please, and determine the corresponding values of 
X. Of course we shall merely obtain the same curve over again. 

It is easy to see from either form of the equation that the 
curve lies entirely above the axis of a: and is symmetrical about 
the axis of y. For from the equation y=x^ we see that y is 
essentially positive, and that two values of ^, which differ only in 
sign, give the same value of y. Again, from x=- + ^y we see 
that y, in order to have a real square root must be positive, and 
that for every positive value of y there are two of x, differing only 



The cm-ve is that given in Pig. (5). 
ment is the same for x and 
y, and the same scale is 
employed in Fig. (3), the 
ciii-ve in Fig. (5) is the 
same as the curve in Fig. 
(3) turned thi-ough a right 
angle about the point 0. 
In fact, the equations of 
the two CHi-ves being y*=iF 
and y=x^, one is situated 
with regard to the axes of 
X and y in the same way 
that the other is situated 
with regard to the axes of y and m. 



If the scale 




§ 1. In arithmetic the numerical difference between two 

quantities is the amount bj which the greater exceeds the less. 
Thus the difference between 5 and 7 is 2 and the greater number 
mity be obtained by adding the difference to the smaller. In 
algebra the amount by which ft differs from a is defined to be b—a 
and this quantity may of course be either positive or negative. 
When we speak of the difference between any two quantities a 
and b, we shall mean the quantity b—a, or the quantity which 
must be added to the former to get the latter. 

^ 2. Suppose we have a series of quantities proceeding accord- 
ing to some given law, e.ff., the fourth powers of the natural 
numbers. Bv subtraction we may find the difference between 
any one and the next following; e.g., the difference between 3* 
and 4* is 256—81 = 175, and 4' may be obtained by adding 17S 
to 3*. These quantities and their differences may be arranged in 



a table thus :- 



1 



1 

16 

8) 

256 



1,2 



1,105 



50 
110 
194 
302 

434 



60 



108 



t7 2,401 
Here the first column gives the natural numbers, the second their 
fourth powers, and the third the differences of the fourth powei-s. 
In the same way that the differences of the fourth powers are 
found, may also be found the differences of these differences. 
These last are written in the fourth column and are called the 
second differences of the quantities in the second column. In the 
same way may be found for this or any other series of quantities 
the differences of the second differences called the third differences, 
and BO on, the process of differencing being capable of repetition 
without limit. It may of course happen that after a certain point 
nil the differences vanish. In the illustration given above it will 
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be found that the fourth differences are all equal, and therefore the 
fifth and higher differences are zero.* It is obvious that the inth 
differences of the nth differences of a series of quantities are the 
m-hwth differences of the series of quantities. 

§ 3. The attention of the student is drawn to the following 
table or scheme of differences : — 



• • • • 

— 5A u_^ • : ; 

— 4A u_ih . ^^u_sh ., : 

— 3A u^2h . ^^u_3h A^u^^.. 

Wo . A2w_,, A%_2A.. 

Amo .2 A^M.;^ 

LiUh A^Wn 

2A fHh . "" A'^% .3 ' A^uo . . 



4A ^4;^ 



AW3/t 



In the first column are a series of equidistant values of x, in the 
second the corresponding values of a function of x denoted by Ux, 
in the third the first differences of the quantities in the second 
column, which may be found by subtraction thus: — 

Am_5^ = W_4A— W_5fe . . . AM_4ft = W_3;t— M_4^ . . . AWs^ = ^4/^ — Wg^ . . . 

In the fourth column are given the differences of the quantities 
of the third column, which are called the second differences of the 
quantities in the second column, and so on. 

§ 4. If we are given the first of a series of values of a function 
and all the first differences, we can, by addition, form the terms 
of the series. Let for instance itoy ^hy ^h • • • denote the series 
of values of the function, and suppose we have given 

Uo, Awo^ ^Uh, i^thh • • . 

Adding Aw© to Mo> we get w^; then adding Auh to Uh, we get tt2A> 
and so on. 

* We do not say that there are no fifth or higher differences, but that the 
fifth and higher differences vanish or are all zero. This may appear a distinction 
without a difference, but it will enable the student to better appreciate some of 
the formula which occur later on. 



Or, if we have given «<,, Auoj A*«o - - . A"Koi and all the n + 1th 
differences, we can form the series Ug, u/,, Uj/, . . . hy addition. 
For by addiog in aucceasion the «+ 1th differences to A"Wo we form 
the table of nth differences, then adding these in snccession to 
A''~'Wo we form the n— 1th differences, and so on. 

Ex.: Given Wo=l, Aw,= 15, A*M(,= 50, A»Mo=60, A'tt„=24 
fof all values of n, find Ug, u/,, u^ . . . 

Here by repeated addition of 24 to 60 we get the column of third 
differences 60, 84, 108 . . .; then adding these in succession to 
50 we obtain the column of second differences 110, 194, 302 . , ., 
and so on, finally obtaining for uj,, U3J,, u^ . . . the values. 16, 81, 
256 .... The figures are shown in the table p. 

With respect to the series of quantities Mqj «Aj "sIj Ws* - - -j 
Uo, Auo, A^Uu . . . are called the initial term and initial differences, 
or the leading term and leading diffei'cnces. The leading term 
and leading differences for the quantities u^,, m.^-j:i.a, w,7+2.* ■ . . 
are of course u„a, Au„a> A^a .... 

§ 5. In sections (1) to (4) we have regarded the differences of 
any function as found by subtraction from a series of equidistant 
values of it supposed given. It will he useful to consider them 
from a somewhat different point of view. 

Let us take any function of x denoted by «, and find the 
increment of the function corrtsponding to an inci-ease h in the 

iahle x, in other words, find the difference between m.^ and 
tix+h- This is obviously u^+k — «*■ The expression Wa+s— 1*^ is 
called the first difference of u^ corresponding to the increment A of 
the independent variable x and is denoted by Awj;. In tlie same 
way Aux being a function of a^ we may find its difference correspond- 
ing to the increment A of x. This is called the second di&crcnce of 
the original function and is denoted by A^u.r- Similarly, the 
third, fourth, &c. differences may be found, denoted by A'm,, 
A V 

Ex. : Take the function a* ; then 

Ai!*=(x + /l)* — 3.'* 

= 4^.^h + G.vVi* + -Ixh^ + A* 
AV = 4(j-|-A)V,-|-6(,* + /0Vr + 4(^ + A)/.HA^ 

= 12jM* + 24j'^^-|-14A', 



14 
Similarly A»^ = 24a?A» + 86A* 

AV=0. 
If h is taken = 1^ we have 

82a^=12^4-24r4-14 

S3^=24r4-36 

S*^*=24. 

from which, putting a? = 1, 2, 3 ... we form the table of differences 
given, p. 

It is to be noted that in forming the successive differences of 
any function, the increment of x, which might also be called the 
difference of x, denoted above by A, is taken to be the same 
throughout. In practice it usually has the value unity, but as 
the investigation of the more elementary theorems of the subject 
is not simplified by assuming A=l, we shall, unless it is otherwise 
stated, suppose its value unrestricted. 

§ 6. The following illustrations of the subject are worthy of 
attention. The demonstrations follow directly from the definition 
of a difference, and they should be tried as exercises before reading 
the proofs. 

(1). If a is a constant quantity, that is a quantity independent 
of a?, 

(2). If Ux and Vx are two functions of x, 

(a) i^[Ux±Vx)^^Ux±i!^Vx 

(y8) ^UxVx = WajAVa- + Vxli^Ux 4" Aw^ AVa? 

Vx VxVx+h 



(8). Alog^ = log(l4-J) 

(4). A*a*=(a^-l)»a* 

(5). A^x^ = \nh^, 
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(1). Aa«ia.=att.^+fc— flftti 

= aAux 

(2). (a) ii{u,r,±V:r)=:U:e^k±Vx+h'~(Ux±Va) 

= Attar ± AV;g 
(fi) ^Ua,Va> = tta?+ Al?« + A — ^arV* 

= MajAVa. 4- t?arAttar + Atta?Al7x 



(7) 






_ t?x(Wa? + Atto?) — Wx(Var + AVa.) 




^V»+fc 




_ VafAttj. — ttj.Al?a? 




t^xt^ar+A 


(8). 


Alog ^= log {x + A) —log X 




='»«'4-* 




=i.g(i+^) 


(4) 


Aa*=a*+*— a* 




= (a»-l)a» 




A»a»=A{a»-l)fl'' 




= (a»-l)Afl« 




= {a»-l)(a»-l)o» 




= (a»-l)»a« 


Similarly 


A»o«=(a»-l)»Afl» 




= (a*-l)»o. 


and ao on. 





(5). This ia a particular case of the following more general 
theorem. 

Let ttjp=A4;f*+Ba?*"' + Ca?*~*+ . . ., a rational integral alge- 
braic function of x of degree n» Tben^ if r<n^ A^^** is a rational 
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integral algebraic function of x of degree n— r in which the term 

involving a?'*"'* is An.n— l.n— 2 . . . w— r + l<r'*~**A*'j if r-=^ny 
A*'aj'* = AjwA'* j if r>n, A**aj'*=0. The proof of this theorem is as 

follows : — 

^-Aa?**— Ba?**"*— . . . 

=Aw^-iA4-Biar'*-2+Cia?'»-3+ . . . 
Differencing again we obtain in the same way. 



similarly, 

b^ux = Aw . n^ . n^2a?»*- W + Bga^"-'* + Ca^;'*-* 4- . . . 
and so on. Thus if r<n we obtain on differencing r times 



A**Ma?=An.n— 1 . . . w— r + la?'*~*'A**-}-B^'*~*-'*4- . . . 
Again, differencing n times we get 



A'*Ma.=Aw.n— l.n— 3 . . . l.A**; 

and this quantity being independent of Xy all higher differences 
vanish. 

§ 7. To express A*^;,. in terms of Wo? and its successive values 

'^x-'t-h} **a?+aft .... 

From the above examples it might be inferred that we should have 
generally 

where pi, p^ » » - are the numerical values of the coefficients in the 
expansion of (a— a?)**. Let us assume, then, that this law holds 
for n, and examine whether it will hold for w + L .Differencing 
we obtain 



Now, if we multiply both sides of the equation 

{a—x)'*=a''--pia''~'.T+pia"~^x^—p3a"~^a?+ . . . 
by a— a:, we obtain 

Thus the coefficients in the aevies for ^"'^'ux are the same as 
those in the aeries for (a—x)"*'. If, then, the law holds for n it 
holds for «+l. But we have seen that the law holds for Ak^, 
therefore it holds for A-Wj-, and therefoi-e for A^Mj., and so on 
universally. That is, we have 



<*+nfe — n«*+S=I/i+ - 



-1 



-1.1 



§ 8. To express «»+ma in terms of Uj- and its first n lending 
differences, i.e., in tei-ras of «», Am^, A^Uj. . . . A"Mj. 

= u, + Aa,+ A(u^+A«,} 

= «,+3A«.r-i-AV- 
Similarly u^+s/i = «j.> + 3Amj, + SA'Uj, + A^Ui, 

and in the same way as in the preceding proposition it may be 



I 



thown that 



. + i,Au.+ 



-1, 



i.a 



, +»A-'«,+a"«, 



(3) 



§ 9. In § 2 wc have supposed a series of quantities ^fHfeerfjV;^ 
accurding to sume ffireii /an: We may, however, take a series of 
numbers at random and obtain from them their first, second . . . 
differences by subtraction. These quantities and their differences 
may still be symbolised as in the table of § 3, the suffix ith of any 
one of the quantities a„it merely denoting its position in the series. 
In sections 4, 5, 7 and 8, no assumption was made as to the nature 
of the law of formation of the quantities considered, and they 
therefore apply to a series of numbers chosen at random and their 
differences. Of such a series of numbers it is perhaps prefei'able 
to say, not that they are numbers chosen at random, but that the 
law of their formation is that of arbitrai-y selection. 



J 



CHAPTER II. 



Interpola 



§ 1. It frequently arises, in the calculation of tables and in 
physical and statistical researches, that certain I'alues of a function 
are found by cxperimentj observation, or calculation, and it is 
required to obtain from them approximately other values of the 
function. For example, we might be given as the result of 
observation the unmber of persons attaining the ages 15, 20, 
25 . . . out of 100,000 alive at age 10, and require to find 
approximately the numbers attaining the intermediate ages 11, 
12, 13, U, 16, 17 

To denote this process of approximating to unknown values of 
a function or terms of a series by means of other values which are 
given we make use of the expression interpulation. 

When we do not know the general form of a function w^ — i.e., 
either the algebraic expi'ession for it or the curve which represents 
it- — but are given merely a series of values of it corresponding to 
a series of difi'erent values of a: — that is, are given merely a series 
of isolated points on the curve which represents it — the problem 
of finding the values of the function for other values of a: or of 
determining other points on the cui've is clearly indeterminate, 
since between any two consecutive given points we may draw any 
number of lines we please, and thus we may get an infinite 
number of curves passing through the given series of points. It 
is therefore necessary to make some assumption in order to have 
a problem with a definite solution. This assumption is usually 
made as follows :— 

§ 2. In the case of most of the tabulated functions — e.g., 
logarithms — the successive orders of differences rapidly diminish, 
and this will usually be found to be the case with the functions 
with which the actuary and tlip statistician are concerned. 

Wheii m particular values of a function are given, and it is 
required to find other values or a general expression for the 
function, it is usual to assume that the function can be completely 
represented by a rational integral algebraic expression in .t of degree 



m—\. This, as we have seen (Ex. 5, p. 15), is equivalent to 
assuming that the with and all higher differences vanish. Pro- 
ceeding on this hypothesis we shall first consider and illustrate 
two particular problems in interpolation, and afterwards (Chapter 
rV.) give a few more general theorems on the subject. 

§ 3, I. — Given m equidistant values Moj wji Wa/, ■ ■ ■ u^;izrih of a 
function, to find a general expression for it — that is, to find its 
value for any other value of the independent variable. 

By the equation (2), p. , we have 



fnAKo^ 



A^-f 



n-l.n-2 



1.2 """ ' 1.2.3 

or, using the symbol x to denote nk, 



A=Ko+ ■ ■ 



where it must be remembered that the only restriction on x is 
that it must be a multiple of A. The last term of the series ii 



since under our hypothe 
vanish. 



i A^Mo and all the higher differences 



Equation (3) gives us a general espression for u^ when 't is 

any positive integer. But the form of the expression for «^ is, 
of course, independent of this restriction as to the value of x. 
Having found the general expreesion for m^ when x has any one 
of the values 0, A, 2h, Sh . . . which are infinite in number, we 
infer that this is the expression for u, when the value of x is 
unrestricted. We therefore have whatever may be the value of x, 



= Mo+ r AmoH- 



^■■(^-/') 



A»u.+ . 



.{x-h)...(x-r^2h) 



(S) 



It will be observed that the above expression is of degi-ee m — 1 in 
X, and giving to a: in succession the values 0, k,Zh . . . m—lh, it 
takes the values Uo, ui,, u^ , , . u^-vi- 

The student should bear in mind that this expressio n for Mj 
is exact only when u^, is known to be an expression of the m— 1th 

B 3 



lately, since A"'tio) A'""'''«|, . . . 
ill usually be convenient to 
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degree in at. If, as will generally be the case,' we know nothing 

of the fuuetioii except the m given values of it, we can only 
that the series represents it approxin 
have all been neglected. 

In numerical applications it w 
write the equation in the form 

^, K«-')„ 1(1-0 ■■■0-'"+'), , 

+ ^AMn+ ^ g ^ A%+ ... + — -;^1 A"'-'Uo . . 

We may regard x as the distance of u^ from the first of the 
given values of the function, and A as the distance between any 
two consecutive given values. 

§ 4. Ex. : The Northampton 3 per-cent annuities for ages 31, 
25, 29, 33, and 37, are respectively 184708, 178144, 17-1070, 
16-3432, and 15-5154. Find the annuity for age 30. 

Taking 15'5154 as the initial value of the function, ive have 

A=4, s:=7, and , =1-75. By subtraction we obtain the follow- 
A ■' 

ing table of differences : — 

Mb 15-5154 .Q.,„Q 

■" i5tSn ■'«38--°5iS -0076 ^,, 
«u 17»144 p-^, — OolO 
«„ 18-4708 °°""' 

Whrace w„=15-5154, Aw,=-8278, A-«o=--064, A'((o = '0076, 
A*Ko= — 0022, aud thu expressiou for H^ takes the value 

1'7!: 



15'515-l, + ]-7oX'8378- 



< ■06-1 



1-75 X ■75 X -35 



X 0076 - 



l-7i 



;^75y85xl25 
24 



■0022 



= 16^9216. 

(The value given iu Jmies un Annuities, voi. i, p. 244, is 16'9217.) 

lu the above aulutiou the initial value of Wj was taken to be 
that correspouding to age 37 and increments of x to decrements 
of age. The problem might equally well have been solved, start- 
ing from age 21 as the initial point and taking increments of jr 
to correspond to increments of age. If this had been done, the 
signs of the odd differeuees would have been reversed. We should 
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have had .r = 9, Wo = 184708, Awo=--6564, Ahio^-'O^l, 
A3wo=— -0054, and A^Mo=--0022. 

§ 5. If we choose to denote the initial value of u hy u^, the 
equation (3) becomes 

Uk^x = WA:+ ^ Am* H- —^j^Dl^Uu + . . . 

This is quite obvious when we bear in mind the remark made at 
the end of § 3 as to the meaning of x, 

§ 6. II. — Suppose the given values of the function do not 
form a complete series of equidistant values in consequence of the 
absence of one or more terms. Let it be required to supply these 
missing terms. 

Let Uq, Uk, iJhJi . .. . denote the complete series, of which m 
terms are supposed given. If one term only is wanting, the 
number of terms in the complete series is m 4- 1, and, to find the 
missing term we put A"'wo=0, or, by equation (1) 



In this equation all but one of the m-\-l quantities Uq, u^jV^*" Umh 
are known, and it is therefore a simple equation to find the one 
required. 

If two terms are deficient the number of terms in the complete 
series being m4-2 we must also put A*"w;i=0, and thus get, to 
find the two unknown quantities, the two equations 



. m,m — 1 . , ,x^ ^ 



ua^h—rnumh-h — j-g— ««sr=i;i— . • • +(— 1)**m;i=0; 

and generally, if r terms are deficient, the total number of terms 
being m-\-r, the r equations to find them are got by equating 
separately to zero the series for A,^Uo, A*'*w;^, i^^thth • • •> A*"w;3ifc. 

Ex,: Given Uq, u^, u^, u^, u%y find Ui and U4. 

Putting A*Wo and A*Wi=0, we have 

1*6 — 5m4+ 101*3—10^2 4- 5t«i — tto=0 
We — 5«*6 4- IOW4 — lOttg 4- 5ti2 — Wi = 
From which we obtain 

«*4 = it{^"" 1 S^ + "^^3 4- 24t*6 — 5 We} 
u, = ^{2wo4- 15wa— 10w3 4^3w5— We} 
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§ 7. The method given in the preceding section will, of course, 
always succeed. In special cases, however, other processes of 
calculation may be employed with advantage. For these no 
general rule can be given. The following is an example {vide 
Journal of the Institute of Actuaries, vol. xv, pp. 894, 895). 
If Mo = 100,000, it; = 97,189, 1*8 = 96,720, ^9= 96,195, find «„ 

Here we assume that 4th differences vanish, and therefore 8rd 
differences are constant. Reverse the series and let the terms in 
reverse order be denoted by Vo, t?i, V2 • • • t^g. By subtraction we 
form the following difference table : — 

vo 96,195 .„. 

V, 96,720 i^Q~56 

V2 97,189 ^^ 

But by equation (2) or (3) 



, . , w.n—l .„ , n.n— l.n— 2.- 
Vn=Vo'\-nAvo -f ^ ^ A% 4- ^ A'vo 

whence, putting n=9, we have 

100,000=96,195 + 9 x 525-86 x 56 + 84 x A%, 

which gives A% = 13*05 . 

Having now obtained the leading quantities i;o= 96,195, Aro=525, 
A%= — 56, A^ya?= 13*05, the table given below is formed by 
addition. 

52500__.^.^ 
46900 _J5.Xr 1305 
42605 _g| J 
39615 _f^.^" 
379-30 %l^ 
375-50" ^12 

^^^•^^ 22-80 
40705 ii^ 

442-40 ^^^^ 



U9=ro 


96,195 


M8 = l'l 


96,720 


«7=»j 


97,189 


«8 = fs 


97,615 


Mb = »4 


98,011 


M« = »S 


98,391 


M3=»« 


98,766 


U2 = Vj 


99,151 


«I = »8 


99,558 


tto=f» 


100,000 
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CHAPTER III. 

Illustrations of the Subject. 

§ 1. To sum the series formed by n consecutive equidistant 
values of a given function of x. Let Ux denote the function, and 
^r> ^+A • • • Wr+Ti^/i the series of terms to be summed. Let Vx 
be such a function of x that Avar=^.r- Then 

Wr + ^r+^H- • • • -\-V'rJf.^\h 

= At7r + Aiv+7»+ . . . 4-Atv+,-^i/, 

So that the sum of the series can be found if the expression for Vx 
can be found. 

Ex. : Find the sum of 



(flr4-6)(arH-l4-6) . . . (ar + w — 1 + 6) 



-\-{ar-\-\-\'b){ar'\-2-\-b) . . . (flr + m + 6) 



4- . . . +(fl.rH-7i— l-|-6)(a/'4-w-f 6) . . . (a.r-l-w + w— 34-6). 



Denoting the function [ax + b){ax -\-\-\-b) ... {ax + m — 1 -f J) 
in which the number of factors is m by wi*"^, and taking the incre- 
ment of X unity we have 

Zuf^=[axTl-^b){ ) . . . (oaj-hw + A) — (aa?4-6)( ). . . (flaj + w— 1 + 6) 



= («..'c+l +6)(fl.aj + 2 + 6) . . . (flaj + w — 1 +6).am 

(m-l) 

Therefore, changing m into m + 1 and a; into aj— 1, we have 
^ hu^]X^=u^r^] that is, — JL^4^1V> is afunction whose 



a(m + l) ^-' ^ ' ;'«(m + l) 

difference is «^$?'*\ The sum of the series is therefore 

a(w+l) ( ) 

= -7-^-^ I (ar + n3l + 6)( ) . . . (fl.r + n + »i--l + 6) 
a(w+l) ( 

— (a.r^+6)( ) . . . (a.r + m— 1+6) 1. 
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§ 2. To find the sum of a series of equidistant values of any 
rational integral algebraic function of x. 

Let the terms of the series to be summed be denoted by 
Woi Wii W2i • • . tin-i} and suppose the increment of x in passing 
from any one to the next is h. Let Sn=«^o + Wi + t/2+ ... +w«-i, 
So being equal to zero. Then, A having reference to the incre- 
ment A in a; or 1 in w, 

ASn=M»; 
Therefore A^Sn=^Un 

A3S^=AX, &c. 
Now, by equation (2), Chap. I, 



S«=So4-wASo4---^A%+ . . . 



. w.n— 1 . . n.n— l.w— 2.„ 
= wt^oH , .rt— AwoH TH A«t^o-h .... 

Substituting in this equation their values for Uq, Auq, A^o • • •> 
we obtain S„ or the required sum. 
Ex. : Find the sum of the series 



K«+(K4-A)3 + (K4-2A)3+ . . . +{K4-n-lA}» 
By the formula just found the sum is equal to 



nK« 4- ^•-^^1 (SK^A 4- SKA^ + A^) 4- "^ * "" -3'^*—^ (6KA« -f GA^) 



. n.n — l.n— 2.71—3^,, 
H = brr 



As a particular case, putting K=A=1, we have 

§ 3. The differences of the various positive integral powers of 
X) in which the increment of x is taken to be unity and x is put 
equal to zero after evaluation, are called the Differences of 
Nothing. The general symbol for them is S'^0'*. To obtain 
S^O^y x^ is differenced m times, and in the result x is put =0. 

These numbers are of importance in the study of Finite 
Differences, and we shall see a use for them later on. We proceed 
to investigate a formula by which they can be easily calculated. 
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By equation (1), Chap. I, putting M^=aj", we have 



7H Vfl ^~ 1 

S"»a;**=(aj + /w)**--m(aj-|-wi — 1)**H — '—- — (aj + m— S)**— ... (a) 

I? 






=inS^""* !»*"*, (by equation {a)). 
But since we have always t^a?+fe=Wa;^- At^j., therefore 

therefore S'^O** = w { S*«- ^0»- ^ + S'^O^" » } . 

By means of this equation the Differences of Nothing may be 
tabulated without difficulty, as in the following diagram : — 



••.. 


01 


02 


0» 


tH 


o» 


o« 


s. • 


■■•••1.. 


1 


1 


1 


1 


1 


s. 





-■■■K 


6 


14 


30 


62 


s. 








6 


36 


150 


540 


84 








- 


•-•24 


240 


1,560 


85 














.120. 


1,800 


8« 

















'■••..720 



Each first difference is unity, since 1**— 0"=1. Again, if the 
diagonal of the diagram passing through the right bottom comer 
be drawn (indicated by the dotted line), all the terms below it are 
zero, since S^i»'*=0 if m>n. The rows may be completed in 
succession, working from left to right as follows : To find a term 
in the mth row, take the term to the left of it and the term above 
that, add the two together, and multiply the sum by m. Thus 

8202=2(04-1)= 2 
820^=2(2 + 1)= 6 
S20*=2(6 + l) = 14&c. 
Similarly 8303-3(04.2)= 6 

880*=3(6 + 6)=36&c. 
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CHAPTER IV. 
Interpolation FoRMULiE. 

§ 1. Given a series of equidistant values of a function, to 
insert between each two consecutive values w— 1 equidistant 
terms so as to divide each interval into n equal intervals. 

Denoting any pair of consecutive given values by u^ and ujc+n^ 
the terms to be interpolated between them will be denoted by 
^k+i) ^k+2 • • • Uk^n-i' The solution of the problem is obtained 
from equation (3). Taking the form of it given in § 5 of Chapter 
II we have 



n \ATnr 



In this equation putting x in succession equal td 1, 2 . . . w— 1, 
we obtain the values of the terms to be interpolated in the interval 
between uu and wa:+«. 

To calculate each term separately by a direct application of 
this fonnula would be a tedious operation. The process may be 
shortened as follows. 

Let S, 5^, S^ . . . denote first, second . . . diflferences when the 
increment of x is unity. If Zujc^ ^uu . . . are found the terms to 
be interpolated can be formed by addition. 

To find ^U)c we must take the rth difference of u^^j^xy that is 

of the series wjfc + - Awjfc 4- \ » A^^j-h . . ., with respect to «, 



n 



2n^ 

and in the result put aj=0. The only quantities which have to 
be differenced are the various powers of x which occur, and the 
calculations may therefore easily be effected by means of a table 
of differences of nothing such as that given §3, Chapter III. 
Writing the series for u^+x in the form 

X J, . x^^nx\^ . a^— 3naj^4-2n*aj . « 
u,,x^u,+ ~A+^^A^^ ^, A3 

+ J-i ^ + ITS ^ + • • • 



whei-e A, A*, . . . are abreviations for Aw*, A*Ui . . ., and making 
use of our table of the Differences of Nothing we write at once 



S«*=- + ^-;rA> + 






l-Gn+lln^-Gn^.. , l-10?«+35«»-50n3-t-24n\ 



\in* 



f 23n^ . , 



30-140«+210«*-10On> , 






|5n* 



\^n^ 



|5n' 



Suppose for example that n=10 and that six orders of differences 
are retained. Then 

S = -1A-045A'+0285A^-0206625A' 

+ 0161 1675 A= - -01 31620125A* 
ff* = -01 A= - OOgA^ + -007725 A*- ■0066975A'' 

+ 005895225A* 
S^=001A»-00135A'+0OU625A'-O014805A« 
B'= -OOOIA'-OOOISA^-I- ■00O2355A« 
8" = -0000 1 Ai* - ■0000225A" 
S"= -000001 A« 

(compare Journal of the Institute iif Actuaries, vol. xiii, p. 160). 

§ 2. As an example of the problem discussed in the preceding 
section, we might be required to Jind the logarithms to base 10 
of all numbers from 1,000 to 10,000, having given the logarithms 
of the numbers 100 to 1,000, that is the logarithms of the numbers 
1,000, 1,010, 1,020 . . . 10,000. The table of differences A, A' . . . 
would be formed from the given logarithms by subtraction, the 
differences 8, 8", S* . . . would then be found for each interval by 



means of the preceding equations, and these differences being 
found the quantities to be interpolated would be found by addition. 
The quantities ^, A^ . . . change in passing from one interval to 
the next, but the coefficients by which they are multiplied in the 
expressions for S, S^ . . . of course remain unaltered throughout 
the calculations. 

It must not be inferred that this would be the most expeditious 
way of makiufr the interpolations in question. The illuHtration is 
given merely to assist the student in understanding the nature of 
the problem discussed in ^ 1. 

§ 3, Lagrange's Interpalaliim Formula. Given )» values of a 
function which are not equidistant, to find any other value of the 
function. 

Let «D) a* I Hr ■ ■ ■ «A be the given values, corresponding to 
the values a, h, c , . . k of .c. It is required to find the general 
expression for Ux- 

Assuming m^ is a rational integral algebraic function of .n of 
degree m — 1, it may be put in the form 



A(»-i)(«-c) . . . (,»■ 



-*) + B(,,— <,)(*-c).. 
..+K(*-»)(*-4)., 



{'-*) 



where A, B, C . , . are quantities independent of x. For this 
espression is of the m— 1th degree in a; and since it contains in 
independent constants A, B, C . . . K it may be made, by a proper 
choice of these Constanta, to repi-esent any hinction whatever of 
the m— 1th degree in ,v. We have to find what values must be 
given to A, B, C . . . K that the expression may represent u^ . 
Putting ^=fl we see that we must have 

«.=A(a-6)(«-.)...(o-4) 



■■ 


(.-4)(.-o) . . . (»-*) 




Similar values af 


oltainecl fur B, C . . . Tlnu ^ 


cgct 


(x-i)(;t-C 


■ ■ ■ <'-*) 1 „, ('-")("-») ■ 


• («-*) 


•(•-»)(«-' 


. . . („_i) 1 «• (S-„)(S_») . 


. (4-*) 



This is called La<:ra. 



( »-a)(»-i) . 
'»(i-o)(*-i) . 



s I'tti'iiiiilii ti>r luter]Julatiuii. 




29 

§ 4. The following transformation of Lagrange's formula is 
of use in the theory of interpolation. 
We have 

1 4. ^-"^ 1 (^— a).(a7— i) {a:^a){x—b){x—c) 
■*" fl-i ■*" (fl-i).(fl-c)' ■*" (a-i)(a-c)(a-£0 

— ^""^ 1 (^— g)(^— &) , (a;— a)(a?— i)(a?— c) 
"^^■*" (a-6)(a-c) ■*" (fl-6)(a-c)(fl-rf) "*" *•• 

fl— 6( a— c (a— c)(fl— rf) ' ' *i 
and this expression may be shown in the same way to be equal to 

(a— 6)(fl— c)l a—-d (a— rf)(a— ^) * ' ') 

«. fl 11 (a?— 6)(a7— c) . . . (a?— A) 

~ ' ^* y (fl-6)(fl-c) . . . (a-A) 

By means of this equation Lagrange's formula may be written 

_ (-I , ^— ^ , (a?— g)(a?— &) (a:— a)(a?— 6)(a?— c) ) 

"*""M ■*'a"^ + (a-6)(a-c) + («-i)(fl-c)(fl-rf) "*■ •••/ 

+ 

"*• (b^a){b--c) ■*■ (c-fl)(c-i) J 

. !f* ) 

■^ (*-a)(A-i) . ..{*-A)i 

where k denotes the last of the quantities a, b, c . » » and A the 
last but one. 
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§ 5. As a particular case let each of the quantities b, c, d . . , 
exceed that which precedes it by n, i.e., let b^a=c—b= . . . 
szA:— A=n. Then the expression 

+ 7a — T/i. \ + W 



(fl — 6)(fl— c) . . . (A — a)(i— c) 
in which the number of factoi*s in each denominator is r 



n'- (jr l.|r--l "^|^b---2 '" 



(-l)^r r r.(r-l) | 

L- ._ 

Now, the increment of x being n, 

A y _ r.(r— 1) 

iA Ujff — Ux+nr ^*^a:+r— 1.71 H jo Wa;^^__2» • • • 

From which we see that the expression (a) becomes —^ . 

The expression (a) involves the quantities a, b, c . . . sym- 
metrically, and we therefore see that if they can be arranged in 
any order a, fi, y . . . such that )8— a=7— /8= . . . =n, then it 
takes the value A^Ua 

Now take a= 6=n 

c=— n d:=2n 

e=z-'2n f=Sn 



• • • . . » 



Then the expression for u„ becomes 

,x(x'-n*)<a;-2n).. , a^(ar'-n'')(a;«-4n«) . . ^ -a, 
■^^ ^ -'A%-*»+-i p '-A^.^+... 03) 

Again take a=0 6=— n 



• • • • • • 



Then the expresaion becomes 

+ — n^i— ' -^.*K_ia+~ F^s— — -A'm_3„+... 

Adding these two expressions for iix we obtain Stirling's 
Interpolation Formula, viz. : — 



!4n* ^" on' 



3»> 

1A>^„ + «_„) + ... 



§ 6. Stirling's interpolation t'ormula gives in another form a 
solution of the problem considered in § 1. If we assume that 
the mth and higher differences vanish it may be shown that the 
two give identically the same values for Ux. These conditions 
however will not be exactly satisfied in any actual case, (except of 
course where the function may be represented exactly by an 
expression of the m— 1th degree in k). The results obtained 
from the formulse of sections 1 and 6 and from other interpolation 
formulas will therefore generally differ to a greater or less extent, 
and it becomes a question of importance what formula may be 
most advantageously employed with a view to making the calcula- 
tions as short and the results as correct as possible. 

In searching for an interpolation formula to subdivide any 
interval it seems natural to look for one which would involve the 
given values of the function on either side the interval symmetric- 
ally. For instance to interpolate in the interval between Wo and 
Mk we should look for a formula, which would involve u^ in the 
same way as «„, «_„ in the same way as Ka„, and so on. It is 
obvious that Stirling's formula is not strictly appropriate for this 
purpose, On referring to the difference table given p, 12, it will 
be seen that we require an equation of the form 

«,3=AK + «„} + BA«o + CA^Cuo+u_„) + DA«i/_„-h . . . 

The quantities A, B, C, D, E . . . must be functions of x and 
n such that the calculations might equally well he performed in 
reverse order, with the series of given quantities . . . u_j„, w_„, 
Woj "hi "sh . t^sB ■ ■ ■ reversed. Now if we reverse the order of the 
given quantities the signs of the odd differences A, A' . . . will be 
changed but the differences otherwise will remain unaltered. We 
therefore see that when a: is changed into n—x, A, B, C, D, E . , . 



must remaiti numerically unaltered, but the signs of the f 
coefGcients B, D . , . must change. 

A formula of this kind may be found without difficulty by 
the method of § 5. Referring to it and putting 







_(5-»), 



/=- 



. . we get 



»(J-„.)(.,-3„) , .-(J-.-(:r-2n)(»-3.) , 



:4«.' 



;5«' 



Adding this to equation (/3) of § 5, we obtain 






lA"(«.+»_.) + 



.(«-«) a.r- 



A%., 



=!(«.+«.)+ -J— A».+ ^2^, ,-„,._.,, ,3^, 

"With reference t» the iiiteiTal between u^ and m« , the quantities 
Aug, JA'(mo + **_k)j A^«_„, . . . are called "central differences." 
In any interpolation forniiik in whieh m, is expanded in a aeries 
of differences of ascending order, the temi involving the difference 
of any order is called the equation of that difference. For 

instance, in the formula just given, the term i ~ i'^'("o"*""-") 

is called the equation of the second diffeience. 

§ 7. As an illustration of Sections (1) and {6), let us take from 
the Institute of Actuaries' Text-Book, Part I, Interest, p. 167, 
the values of log,„(H-i) for 100/=6tV, ^^, 7^, ■ ■ ■ 10. and 
deduce from them the values of log|o(l+0 for 100i=7-pj, 7-i-|[, 
7H' ■ ■ ■ *iV; *'^^* '^' *^'^ 'M* '^^ interval between 7^ and 8^. 
Retaining 13 tigures in the logarithms, and omitting the 
initial zeros, we obtain by subtraction the following difference 
table :— 

Table A. 



, 383 761 277 363 
314084642 616 
33» 261 204 729 
364 292 656 S67 
3H9180 6G030'1 
413926851 582 



25* 

US 323 365 154 

MS 176 562 213 

25031451638' 

24888004 037 

24 746191 278 



148 524085 
-146802 041 
146110 676 
143 447 501 
141 812 766 



33 

Let us first consider the formula of § 6. Retaining fourth 
differences, it may be written 

+ ^' '^^-»+ »? lA'("-. + «-t.) 

from wtichj by the method of § 1, we obtain 



*= ion S'^'<"'+ "-) - iTSSi'^'''-"- 40000 ■ l"^''"- 



'100' 



n4A'(._..+. 






■iAV. 



whence Bu„ =3 524 33+ 200-559 187 5 

2^= -1466 399-614375 

S'ao= 1703568 25 

S*ut= —2943 5 

The above are the exact values of Bu^, S^o ... as deduced 
from the data. If, instead of these values, we use approximate 
values di, rfi, rfj, df, obtained by omittinp; some of the last figures 
of the decimal parts, and if E|, E^, Ea, E, are the errors thereby 
introduced into 8, S*, S=, S', ao that rf, = « + E,, rf,=S^ + Ej . . ., 
the value found for u^ will be 

^+,,,+ 'fcd>,.+ . . . + ^.l^-i)(^;-2)(..-3) ,^ 

bo that the error thereby introduced into Ux will be 

The coefficients of E,, Eg ... are gi-eatest when x= 10 when 
they have the values 10, 45, 120, 210. If, then, we retain two 
decimal places in Bua, and three in each of the others, the greatest 
error introduced will be less than 10x-005 + (45 + 120 + 210) 
X 0005 = -2376. 

Let us then take the decimal parts of S, £* ... to be -56, -614, 
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respectively. These 
E,= -000375, E3=--0O025, %=-Q005~m that the eiwr intro- 
duced into w,o is -008125-1- •016875--03-|--105 = -l, and the 
maximum en-oi- introduced is less than ■00812a + 'U16875-|--03 
-|--105 = -16. 

We now proceed to the interpolation b. ^Ub is first written 
down (see Table B), and then 8^W[, S^Uj . . . formed from it by 
i-epeated addition of S*u=17'057. These are written on every 
fifth line. The value of Sht^ is then written above S'uo, the two 
added together, and the sum — i.e., S'u, — written above S^Wi, then 
S'ui and S^i are added, and the sum — i.e., S^ua — written above 
S'lia, and so on. In this way the second differences are foimed. 



and then by a similar proces 
quantities Ui, «„ . . . um. 



the first ditFerenccs, and finally the 



S 524 234 200 56 

13 533 700-386 

I 70a'568 

316 608 876 716-56 

2 522 767900-946 

18 535 402-95-I 



326 691167583-625 

2 516 919 698740 

18 54a 183-796 

1 687-853 

»29 20S087 232'365 

2 515 461 882-536 

18 643 871-649 



319 13l644617-50fi 

2 521 303 308-900 

18 537 lOa-579 



321652947 921-406 

3 510 840406-479 

18638799-261 



324172 783327-885 

2 518 379 205-740 

18 640493-000 

1 690'796 



331 723 549 1U'901 

3 514 005 754' 18; 

18 545 556-559 



334237 554 869-086 

2 512 551 310744 

18647238-526 

I 679-034 

336 750 106 179-830 

2 511 098 549-270 

18 648 917550 

1 676081 



If we go back to the formula (4) of § 6, we notice that when 
a:=nthe aeries takes the value «« whatever values wc assign to 
A'(«o+«_n). i^«-» ■ ■ ■ ; and therefore from the fact that Um 
comes out correct, except for the predicted error -1, we may not 
infer tliat the calculations arc free fi-om error, though we may 
infer that the additions of Table B have probably been correctly 
performed. 
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§ 8. The method in which the computations involved in any 

interpolation may most advantageously be effected will depend on 
the particular problem under consideration. Two points may be 
noticed. 

(1). The equations of the differences diminish as the order of 
the differences increases. 

(2). That when, as in § 7, the calculations are effected by 
continued addition of differences, there ia an accumulation of 
error, which frequently increases mth the order of the difference. 
For instance, in § 7, an error in S'zio is magnified 10 times in w,o, 
while an error in S^uo is magnified 310 times in ui„. Owing to 
this it is often necessary to retain in the higher differences several 
figures more than are wanted in the final results. 

For these, among other reasons, it is sometimes profitable, not 
to proceed as in § 7, but to introduce the equations of the higher 
differences into the calculations by a different process. 

Let us now make use of the formula of § (1) to perform the 
interpolations already effected in Table B, first obtaining approsi- 
mate values of the quantities to be interpolated by proceeding as 
in § 7, but retaining only second differences, and then correcting 
the results so obtained by adding the equations of the higher 
differences. 

Retaining only second differences, the formulie of § 1 give ub 
_ Am„ 



8»ao=-01AX 



2 52U86S01-675 



= - 145 110 6-75 
= 1 854 889 3-25 
The Table C. (given below) is formed as follows; — 

8u„=2 524 186 201-675 is first written down at the top of 
the table, then by repeated addition of S^ we form S«i, Skj . . ., 
written on every seventh line. «„ is then written underneath Zu^, 
the two added, and the sum wi-itten underneath Su,, and so on. 

The results so obtained require to be corrected by the equations 
of the third, fourth and fiflh differences, which have been obtained 
from their values 



X If 



x{_x-\0){x-2Q) 

g(10)i 

x{x-\Q) { x-2Q){x-m\ 

^(10)* ** 

j(j-10)(j;-20)(j.— 80)(j- 
[SllO)- 



3 174 
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by direct calculation. These are written in order underneath, and 
each of the final results t^i, t^ . . . is then obtained by addition of 
the four numbers above it. We notice that A*Wo cannot be 
obtained exactly from our data, but on the assumption that sixth 
differences are constant, its value is 634. It has been put =640, 
which is suflSciently exact for our purpose. 



Table G, 



2 524 186 201-675 
Uo 314 084 642 516 

2 522 735 094*925 
316 608 828 717*675 
47400*46 
587-48 
10*3 
Ui 316 608 876 715-9 

2 521 283 988*175 
319 131 563 812*600 
79 832*35 
955*31 
16*4 
U2 319 131 644 616*7 

2 519 832 881*425 
321 652 847 800*775 
98 958-85 
1 141*90 
19 
Uz 321 652 947 920*5 

2 518 381 774*675 
324 172 680 682*200 
106 443-14 
1 182*77 
19*2 
Ui 324 172 788 327*3 

2 516 930 667*925 
326 691 062 456*875 
103 948*38 
1 110*62 
17*5 
U^ 326 691 167 533*4 



2 515 479 561*175 
329 207 993 124*800 
93 137-74 
955*32 
14*6 
U^ 329 208 087 232*5 

2 514 028 454*425 
331 723 472 685*975 
75 674*42 
743*85 
11 
Uj 331 723 549 115*2 

2 512 577 347*675 
334 237 501 140*400 
53 221*57 
500-40 
7*2 
U^ 334 237 554 869*6 

2 511 126 240*925 
336 750 078 488*075 
27 442*37 
246-29 
3*4 
Mg 336 750 106 1801 

2 509 675 134175 
Mio 339 261204 729*000 



§ 9. When a single interval only has to be filled up, it will 
generally, as in the example we have considered, be far easier to 
calculate all the leading differences, and then form the table by 
addition, as in § 7. But if a large series of intei-vals have to be 



filled np, the pmceases may woinetiraes he. shortened either by 
computing the equations of the higher differences directly or 
by the use of various artifices to introduce them into the calcula- 
tions, anbaidiary tables being formed to shorten the work. {Vide 
paper on Interpolation by Mr. Woolhouse, Jourjtal of the Institute 
of Actuaries, vol. xi, p. 61.) 

The figures of Table C afford a useful illustration of the 
gradual closing of the approximations as the successive differences, 
are taken account of, and it is for this reason that the table has 
been given in the shape in which it stands, 

§ 10. On comparing the iatei'polated quantities found in 
§§ 7 and 8 with the values of the quantities as given on p. 167 
of the Institute Text-Book, it wiU be seen that the i-esults of each 
interpolation are correct, or nearly so, to the last figure retained, 
the greatest error being about "5 in this figure. If we calculate 
for Table C the equations of the 6th differences, we shall fiud 
the greatest of them is about '7. For the formula used in § 7, 
the greatest value of the equation of the 5th difference is about 
"6. We therefore see that the results obtained by the use of 
formula (4), retaining only 4th differences, are at least as accurate 
as those given by formula (3) retaining 5th differences. It ia 
easily seen that the scries (4) generally converges more rapidly 
than the series (3). 

* § 11. It is not to be supposed that all interpolation formulae 
are founded on the assumption of § 2, Chap. II. Mr. Sprague, 
in a very interesting paper (Journal of ike Institute of Actuaries, 
vol. sxii, p. 370} has given formulse for interpolation framed on a 
different basis, which deserve attentive study. The hypothesis 
on which he proceeds may be explained as follows : — 

Let the ordi- 
nates Pp, Qq . . . 
represent given 
values yn, y, . . . 
of a function be- 
tween which it is 
desired to interpo- 
late other values. 
" The problem ov 
"interpolating 
" between y^ and 
" ys ia the same 
" as that ov drawing a curvd line between the points fl and S, 




J 



se 

" and in order to get a satisfactory interpolation, it is necessary 
" that this partial curv ahoud join on smoothly to the ajacent 
" partial curvs, namelyj Q R on the one side and iS T on the 
" other side." 

In the formulje for interpolation hitherto considered, this 
smoothness of junction is obtained, with more or less success, by 
drawing the partial curves so that any one, as R S, forms a 
portion of a cui-vc through R, S, and adjacent given points. 
This partial curve approximates to a portion of the parabolic curve 
of degree ni — 1 passing through all the given points, supposed rn 
in number, coinciding exactly with this curve when m — 1 orders 
of differences are retained, 

Mr. Sprague secures this smoothness of junction by drawing 
the partial curves so that any two adjaeent ones, as Q R, R S, 
have, at their common point R, contact of the second order with 
the quartie parabola passing through R and the two adjaeent given 
points on each side, P Q and S T. Each partial cui-ve has there- 
fore to satisfy three eonditiona at each extremity, six in all. 

Let us find the equation to the curve R S, on the assumption 
that the given ordinates are equidistant. Taking r as origin, it 
will be of the form 

y = I/a + n,i; + iiE= + C^ + rf3r' + Ar= I 

Now let us form the table of differences for the quantities 
Voi y>j Vi ■ ■ ■ Vi- Denote the differences of y^ by A, A", A' . . . 
and those of j/i by Ai, A[^, Ai' . . . 

Referred to p as origin, the equation of the quartie parabola 
through P, Q, R, S, T, is 



^ ^(»-l)(»-2)(a;-3) ^. 



Transferi-ing the origin to r, this becomes 

(»+a)(»+i)»(»-i) ^. 
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The equation I, referred to s as origin, is 

+ ^{6 + 3c + 6rf4-10e}+ ... . Ill 

and, by observing the form of Equation II, we see that, referred 
to s as origin, the equation to the quartic parabola through 
QRSTUis 

y=y3+.^(A,-ffAiHiA,«-TVAi^)+^(iA,HiA,«-3jVV)-f..IV 

Since I and II have contact of the second order at R, we 
obtain, equating constant terms and the coeflScients of x and ^, 

«=A+|A*+iA«— jVA^ (a) 

6=|A''+iA»-^A^ (^) 

Similarly, since III and IV have contact of the second order 
at S, 

a+J + c+rf+c=y3-ys=A + 2A*-'+A» .... (7) 

a + 2« + 3c+4</+5c=A,+fA,''+|A,'-TVV 

=A+fA«+VA='+iA^-TVA» . . (S) 
6 + 8c+6«i+10e=|A,''+|A,8-^A* 

^^A^+A'+^A^-^A' . . . . (0 

he five equations (a), ()8), (7), (S), (t), which are the same as 
those obtained by Mr. Sprague by a somewhat different process, 
determine a, i, c, rf, e, in terms of A, A^ . . . A*. Solving them, 
and substituting in Equation I, we get the required interpolation 
formula. 
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EXEECISES AND EXAMPLES. 

1. Show, by actual addition, the succeeding terms of a series 
the first of which is u and the successive differences of whose 
terms are Am, A%, A^u .... 

** Institute of Actuaries* Exam.", vide Journal iii, p. 274. 

2. Point out the analogies which lead us to infer that for 
the nth term we may write tt»=t*-f wAw H — -f-o" A^+ . • . and 

for the nth difference A'*t*=w»— wti^-i -f — , q- «*»-2— • • •• 

" Institute of Actuaries' Exam.'', vide Journal iii, p. 274. 

3. Find 8**^* when x is variable, the increment of x being 

unity. 

** Institute of Actuaries' Exam.", vide Journal xviii, p. 378. 

4. If Wa? be any function of a? of n dimensions, prove that 
A^Ux is constant ; and hence show how to form a table of cubes 
of natural numbers expeditiously. 

" Institute of Actuaries* Exam.", vide Journal xxii, p. 65. 

6. Investigate the expressions : 

(a) for Ua,+n in terms of Aw^j. 

()8) for A^Ux in terms of Ux and its successive values. 
" Institute of Actuaries' Exam.", vide Journal xxii, p. 65. 

6. Investigate an expression for A^Ux in terms of Ux and its 
successive values. 

Using the formula thus found, if in the series 1, 6, 21, 56, /c, 
252, 462, &c., the sixth differences vanish, find k and sum the 
series to 10 terms. 

" Institute of Actuaries* Exam.", vide Journal xxi, p. 224. 



41 

7. Given t^a,= 100,000, 8i=-490, 83= 93, 83= -25, 84=8, 
§6=0, find the first eight teni[is of the series Uof, Wa?+i, &c. 

" Institute of Actuaries' Exam.'*, vide Journal xxii, p. 65. 

8. Find m^ and also u^ when «^5=55, U9=126, ^7=259, 
t«g=484, Ug=S37, and A^ is constant. 

" Institute of Actuaries* Exam.*', 1881. 

9. Given the following values, construct logs 7'1, 7*2 . . . 7*9, 
8, to seven places, and explain why the last diflfiers from log 8 as 
obtained from an ordinary table of logs. 

log 7= -8450980 
Si= 61603 

82= - 861 

83= 24 

84=- 1 

" Institute of Actuaries' Exam.", yide Journal xzi, p. 224. 

10. Given log235=2-8710679 

log 236=2-3729120 
log 237 =2-3747483 
log 238=2-3765770 
Find log 23563. 

" Institute of Actuaries' Exam.", vide Journal xviii, p. 378. 

11. Find log 512, given that 

log 510=2-70757018 
log 51 1=2-70842090 
log 513=2-71011737 
log 514=2-71096812 

" Institute of Actuaries' Exam.", vide Journal xx, p. 139. 

12. Havmg given 

XlOl =2-0043214 
Xl01-5= 2-0064660 
\102 =2-0086002 
\102-5= 2-0107289 
\103 =2-0128372 
X104 =2-0170883 

and assuming that sixth differences vanish find X103*5. 

" Institute of Actuaries' Exam.", vide Journal xxii, p. 65. 
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13. Having given the values of annuities at the following 
rates of interest, namely, at 

3 =15-863 
3^ = 14-941 

4 =14105 
4^ = 13-343 

5 =12-648 

Find the value at 4*328 per-eent. 

** Institute of Actuaries* Exam.", vide Journal xxii, p. 300. 

14. The H^ premium at age 40 is at 3 per-eent '025891 

•024654 
•023517 
•022470 
•021509 
•019811 

Interpolate the corresponding premium at 5^ per-eent, — (1) using 
two of these values, (2) using /ot^r, and (3) using six. 

"Institute of Actuaries* Exam.", 1881. 

15. Show that the sum of the series 

1 1 ._!_ . ^-^ 

m{m^r) ^ (m + r)"(m + 2r) ^ (m + 2r)(m + 3r) ■+■••• ^^ ''V^'*- 

is equal to — . 
^ mr 

"Institute of Actuaries* Exam.**, vide Journal iii, p. 274. 

16. Given every nth term of a series of values, i.e., w^?^ ^x-^ny 
Ux+2ny &c., show at length how the intermediate terms, Wa?+u 
Ux+2 . . . may be obtained by interpolation. 

Given that in the series %., u^+i, u^^^ • • • 

Wo; =9936675-4 
Si =+12767-62 
^2= -3013-725 

83= +422-8247 
84= -34-72847 
85= +1-254221 

you are required to construct the series as far as the term u^r+io* 
What assumption is necessary? 

•* Institute of Actuaries' Exam.", 1876, vide Journal xx, p. 189, 
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17. If Ux be a function of x of the form biX + bf^+Scc. ad 
infin.y show that it can also be expressed in the form 



1_^ ' (l-a?)^ ' (1-^) 
We have 

ft»,+r+i=fti+(m + r)AJi+ ... + ^ '-^ i ^ ^ -' A'»i,+ ... 

so that the coefficient of A*»ii in the term bm^-r^i^'^^^^^ is 



(m + r)(m + r— 1) . . . (r-fl) 



TW 



^t»+r+i 



= ^«M (m + l)(m + 2)...(>»+r) ^ 

Ir 



where P,. is the (r4-l)th term in the expansion of (1— 5?)"^'^+^^ 
Thus the coefficient of A'^fti in u^ is ( vZT ) 



18. Show that 

and hence determine a series of such a nature that the terms 
after the first shall be respectively double the first terms of the 
successive orders of diflference, (t«2=3A*Wi, t«3=2A^o and so on). 

" Institute of Actuaries* Exam.", 1881. 

19. Prove by means of Stirling's interpolation formula that, 
if third and higher differences are neglected, 

__ x(n-]rx) n^^a^ x{n—x) 

Journal of the Institute of Actuaries, xv, p. 392, 

20. If t<4=85, t<5=156, M«=259, W7=400, ti8=585, and 
fourth differences are constant, find 

(i) The differences of 1*4. 
(ii) The value of Wg. 
(iii) The general expression for %, 
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21. Given n„=u^=0, 7(,= 112, ■»,„=9,100, Hso=156,400, 
and that fourth differences ai-e cunstaot, find by means of 
Lagrange's formula the general expression for u^. 

Result, Uj,=x*—9x^. 



22, Let u,, 11^. . . denote a series of quantities, and S, denote 
the sum of the first x of them, So being =0. Having given the 
values of S„j Sj,, . . . S^„, show how to find the values of Tt,, %, 



1(3 . 



- "n.. 



Let 5 denote the difference symbol for the increment unity in 
the suffix, so that S!tj;=Wn.i— m-i, SSj=Sj:+, — Sb=m,+], Let 
A be the difference symbol for the increment n in the suffix, so 
that Anj:=Mjr+„— Mj., AS3:=Sj,^.„— Sj:. Form by subtraction the 
difference table for the series of given quantities Sj, §„ . . . 8„,, 
and let A( = Sa), A^ . . . A"" denote the leading differences. By 
the method of § 1, Chap. IV, we can find BS^, S^So, 8»So . . . S-'So 

in terms of A, A^ But 2S^ = ?(j;+i, .-. ^Sx=Bit^^,, &c., 

so that we obtain the values of Mi, 8w,, S'u, . . . S''"'^,, and then 
by addition of differences we can form the table wi, w^ , . , m™. 

As an illustration take that afforded by Mr. Berridge's gradua- 
tion of the Peerage Mortality Table (vide Journal of the Institute 
of Actuaries, yo\. xii, pp. S20, 221). Here 7i = 10, and the 
difference table for the quantities So, S„, Sa„ . . . is that given 
at the top of p. 231 of the Journal. A = S„ = 99,616,210, 
A"=-18,425, A==-9,898, A*=-186,09e, &c. By the 
equations of § 1, Chap. IV, we have 

t(, = -lA— ■045AH-0285A3— . . . 
8% = -0IA^--0O9A»+ . . . 
S=«, = -001i*- . . . 
&c. &c. 

The first term, w,, ia obtained by Mr. Berridge as follows : — 



S.=S, + .i8S,+ - 



1.2 



-S'S,- 






-\{^-{% 



Thia equation gieea «, when Swi, S'li, ... are found. 
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As another illustration, find Ui, u^y i^ . . ., having given 
S5=l,365, Sio=5,155, Si5=13,370, Sao =28,635, third differences 
being assumed constant. It will be found that t^i = 154, Swi=49, 
S«w, = 10, 33^^1 = 1. 

23. If t^a.= (a?— a)(^— J)(a?— c) . . . to n factors, prove that 
St^^=Si + S2+ . . . +S„, where, for all values of r from r=l to 
r=7i— 1, S,. denotes the sum of the different products that can be 
formed from the factors ^— a, a?— i . . . taken n—r together, and 
S»= 1. (The increment of x is supposed unity.) 

24. If ^(^-1) >■> (^-n + l) .^ ^^^^^^^ ^^ j,^^^^^^ ^^^^ 
that, the increment of x being unity, 

26. The increment of x being unity, we have 

t^a?=^o+^S^o+ "So — ^^2^0+ • • • 
and therefore, taking h'^Ux as the function of Xy 

Verify this by differencing the expression 

UQ'\'xhuQ + -^7, — - S^^o+ • • • 
and making use of Example 24. 

26. Having given t^o=fl^^^-i = *^^-2=^^ S^i^„3=rf^S^^-4=^^ 
prove that 

^ ^ . ^(^+1)' . x {x-\-\){x-\'2) j . ^(^ + l)(^ + 2)(^ + 3) ^ 
^a;=a + a?ft+ ^ ' c + -^ jl^^ ^rf + -^ 1^ ^^ 

fourth differences being constant and the increment of x unity. 

27. The increment of x being unity, prove by induction or 
otherwise that 

Mx=tio + ^OW-l+ -^-o ^ 0*W^2+ — ^ T^ 0't^_3+ 
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Assume this theorem holds for a particular value of x. Increasing 
the suffixes by unity, our hypothesis gives us 



t^ar+ 



l = t^l + ^S2^o+ -^o — ^^^-^1+ • • • • 

If 



Now, 'Wi=t^+Sw_i + S*ul2+ • • • 

whence the coefficient of ^^u^r in the expression for w^+i is 
l+a7H — ~ — - + :^ + . . . H — ^^ '- ^ '- 



and it may easily be proved by induction that this 



r 

L_ 



•X-NOTB. — ^The above theorem follows at once by the method of separation of 
symbols. Its symbolical expression is (1 + 5)*Mo = ( 1 — .j — ^ t«o. 

¥c 28. Prove Briggs* Interpolation Equations (vide Journal of 
the Institute of Actuaries , vol. xiv, p. 79). 

Let S denote the difference symbol when the increment of x is 
unity, A when the increment of a? is 5. Also let 1 + S be denoted 
by E. It may be easily verified that 

^=SE2+S3E + -2S«t 

5 

whence (tT^* = (^^' + ^'^ + '^^^^^ • 

t This is a particular case of the theorem 

+ . . . +«a(i+a)«, 

n being an odd number =2m + l. This theorem was communicated to me by 
Mr. W. L. MoUison^ of Clare Coll., Camb. It may be obtained by equating the 
coefficients of or^ in the equation 

log(l-|?a?)+log(l-g'a?)«log{l-(i? + g)a:+|?g'a:a}, 
and then putting |?«l+5, g'=— 1« 
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In this equation, putting n in succession equal to 1, 2, 3 . . . 20, 
we form the table given [Journal of the Institute of ActuarieSy 
xiv, p. 79) — e.g.y putting n=13, we get 

^^Y'=:SwE26 4.i3gi5E25 4.80-6S»7E24 4.3i7.2Si9E^+ . . . 

whence 

S«i^,^26=^^-{13S»War+25 + 80-6S»7i^^^.24 + 317-2Si»u^+23+...}. 

Writing {SE2+S3E + -28*}»* in the form 8*^*{a2+a + -2}«, where 

a= ^^ we see that it may be expanded in a series of terms, in 

each of which the index of 8 is greater by 2 and the index of E 
less by 1 than in the preceding term, the first term being h'^W^. 

Now, let two tables of dififerences be formed, as at p. 12, one 
of the S differences, the other of the A differences, and imagine 
the second superposed on the first, the two tables being so formed 
that corresponding values of Ux may coincide. Since the interval 
between two consecutive terms of the second table is 5 times as 
great as between two of the first table, it follows that ^^^u^ will 
fall in the 5nth row of the first table below u^- But S~E^t<j, 
= S*Ua.+2», and is therefore in the (4w4-w)th row below Ux — 2.^., 
in the same row in which A"i*.r will fall. And, since in the 
the expansion of {SE^ + S^E^- '28^}** the powers of 8 increase by 2 
and those of E diminish by 1, all the terms of {8E^+ S^E + '2,h^}''Ux 
will lie in the same row. 

29. Denoting Uxy, a function of two independent variables, 
X and y, by a?y, prove that, if third and higher differences are 
neglected, 

1 _ o^ 1 _ 

;ry=i{00 + 0l4l0+ll}+^^^^ {11 + 10-01-00} 

+ ?2!:^ {n + oT-ino-6o}+ ^^^ 

-.y{y-\) 1,_ 



-(11 + 10 + 01 +00)} 4 ^^% ^ .^{12 + 02 + 1-1+0-1 

-(u+io+ol+oo)} + ?fi:i .fcl {Ti +00-61- 10}. 
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30. Find -pr^ {mantissa of logio 41*6}, having given the 
following table of values of the function - {mantissa of logw x}. 



X 

II 


y-40 


41 


42 


43 


40 


• • t 


1,468,439 


1,433,476 


• • • 


41 


1,531,960 


1,494,595 


1,459,009 


1,425,079 


42 


1,558,123 


1,520,120 


1,483,927 


1,449,417 


43 


• « • 


1,545,045 


1,508,258 


• • • 



31. From the following table of endowment assurance pre- 
miums^ 





50 


56 

£2 8 9 


60 


65 


15 


• • • 


£2 4 6 


• • • 


20 


£3 6 


2 17 6 


2 11 8 


£2 7 8 


25 


3 18 4 


3 6 


2 17 11 


2 12 6 


30 


* • • 


3 19 5 


3 7 8 


• • • 



in which the ages at entry are given on the left side, and the ages 
at which the endowments are payable at the top, find the premium 
for age at entry 21, sum assured payable at 58 or previous death. 

32. Prove by means of Chap. Ill, § 2, or otherwise^ that if 
r is a possible integer 



P+2*-+ . . . +n*-= 






where a vanishes when n is infinite. 
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